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Abstract 

The main objective of this paper is to give a rigorous treatment of Wigner's 
and Eisenbud's i?-matrix method for scattering matrices of scattering sys- 
tems consisting of two selfadjoint extensions of the same symmetric oper- 
ator with finite deficiency indices. In the framework of boundary triplets 
and associated Weyl functions an abstract generalization of the i?-matrix 
method is developed and the results are applied to Schrodinger operators 
on the real axis. 

2000 Mathematical Subject Classification: 47A40, 34L25, 81U20 

Keywords: scattering, scattering matrix, i?-matrix, symmetric and selfadjoint 
operators, extension theory, boundary triplet, Weyl function, ordinary differen- 
tial operators 



R-MATRIX 



2 



1 Introduction 

The i?-matrix approach to scattering was originally developed by Kapur and 
Peierls [21 in connection with nuclear reactions. Their ideas were improved by 
Wigner (40j[4T] and Wigner and Eisenbud [42j . where the notion of the i?-matrix 
firstly occurred. A comprehensive overview of the i?-matrix theory in nuclear 
physics can be found in [7J [H] . The key ideas of the i?-matrix theory are rather 
independent from the concrete physical situation. In fact, later the i?-matrix 
method has also found several applications in atomic and molecular physics (see 
e.g. [3 [8]) and recently it was applied to transport problems in semiconductor 
nano-structures [28l ESI [30l EH E21 EH [43l SH [45]. In [26l [27] an attempt was 
made to make the i?-matrix method rigorous for elliptic differential operators, 
see also [3H 133 for Schrodinger operators and [3S1 [37] for an extension to Dirac 
operators. 

The essential idea of the i?-matrix theory is to divide the whole physical 
system into two spatially divided subsystems which are called internal and ex- 
ternal systems, see [40l HTJ [42] - The internal system is usually related to a 
bounded region, while the external system is given on its complement and is, 
therefore, spatially infinite. The goal is to represent the scattering matrix of 
a certain scattering system in terms of eigenvalues and eigenfunctions of an 
operator corresponding to the internal system with suitable chosen selfadjoint 
boundary conditions at the interface between the internal and external system. 
This might seem a little strange at first sight since scattering is rather related 
to the external system than to the internal one. 

It is the main objective of the present paper to make a further step towards 
a rigorous foundation of the i?-matrix method in the framework of abstract 
scattering theory [5], in particular, in the framework of scattering theory for 
open quantum systems developed in [3j [4]. This abstract approach has the 
advantage that any type of operators, in particular, Schrodinger or Dirac oper- 
ators can be treated. We start with the direct orthogonal sum L := A © T of 
two symmetric operators A and T with equal deficiency indices acting in the 
Hilbert spaces Sj and respectively. From a physical point of view the systems 
{A, f)} and {T, &} can be regarded as incomplete internal and external systems, 
respectively. The system {L,£}, £ := Sj © R, is also an incomplete quantum 
system which is completed or closed by choosing a selfadjoint extension of L. 
The operator L admits several selfadjoint extensions in £. In particular, there 
are selfadjoint extensions of the form Lq = Aq © To, where Aq and To are self- 
adjoint extensions of A and T in 9j and R, respectively. Of course, in this case 
the quantum system {Lo,£} decomposes into the closed internal and external 
system {Aq, ft} and {To, M.}, respectively, which do not interact. There are other 
selfadjoint extensions of L in £ which are not of this structure and can be re- 
garded as Hamiltonians of quantum systems which take into account a certain 
interaction of the internal and external systems {A, Sj} and {T, M}. In the fol- 
lowing we choose a special self-adjoint extension L of L introduced in [9] and 
used in [3J, see also Theorem EUJ which gives the right physical Hamiltonian in 
applications. 
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For example, let the internal system {A, Sj} and external system {T, &} be 

,2 

given by the minimal second order differential operators A = + v and 

T = -^5 + V in Sj = L 2 ((xi,x r )) and & = L 2 (R\(xi,x r )), where (xi,x r ) is 
a finite interval and v, V are real potentials. The extension Lq can be chosen 
to be the direct sum of the selfadjoint extensions of A and T corresponding to 
Dirichlet boundary conditions at xi and x r . According to [3, 9J the selfadjoint 
extension L coincides in this case with the usual selfadjoint Schrodinger operator 

7 d? _ jv(x), xe{xi,x r ), 

L = 77 + v. v (x := < 

dx 2 \V(x), x e R\(x h x r ), 

in £ = L 2 (R), cf. Section O 

Let again A and T be symmetric operators with equal deficiency indices in 
Sj and ^, respectively. It will be assumed that the deficiency indices of A and 
T are finite. Then the selfadjoint operator L is a finite rank perturbation in 
resolvent sense of Lq = A n T and therefore {L, L } is a complete scattering 
system, i.e., the wave operators 

W±(L,Lq) := s- lim e 4 * £ e - itLo P QC (L ) 

£ — >±oo 

exist and map onto the absolutely continuous subspace $) ac (L) of L, where 
P ac (Lo) is the orthogonal projection onto F) ac (Lo), cf. |2]. The scattering 
operator 

S:=W+(L,L )*W^(L,L ) 

regarded as an unitary operator in the absolutely continuous subspace Sj ac (Lo) is 
unitarily equivalent to a multiplication operator induced by a family of unitary 
matrices {S'(A)}a6R in a spectral representation of the absolutely continuous 
part of Lq. This multiplication operator {S'(A)}agk is called the scattering 
matrix of the scattering system {L, Lo} and is one of the most important objects 
in mathematical scattering theory. The case that the spectrum cr(Ao) is discrete 
is of particular importance in physical applications, e.g., modeling of quantum 
transport in semiconductors. In this case the scattering matrix of {L,L } is 
given by 

S(X) =1- 2z v / 3m(r(A))(A/(A) + r(A)) "V»m(r(A)), 

where M(-) and r(-) are certain "abstract" Titchmarsh-Weyl functions corre- 
sponding to the internal and external systems, respectively, see Corollary [521 

The i?-matrix {i?(A)},\eR of {L, Lq} is defined as the Cayley transform of 
the scattering matrix {S'(A)}agr, i.e., 

R(X) =i(I-S(X))(I + S(\))-\ 

and the problem in the i?-matrix theory is to represent {-R(A)}agr in terms 
of eigenvalues and eigenfunctions of a suitable chosen closed internal system 
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By the inverse Cayley transform this immediately also yields a repre- 
sentation of the scattering matrix by the same quantities. 

For Schrodinger operators the problem is usually solved by choosing appro- 
priate selfadjoint boundary conditions at the interface between the internal and 
external system, in particular, Neumann boundary conditions. We show that in 
the abstract approach to the i?-matrix theory the problem can be solved within 
the framework of abstract boundary triplets, which allow to characterize all 
selfadjoint extensions of A by abstract boundary conditions, cf. [TO l \TT \ \12 \ fT9] . 
It is one of our main objectives to prove that there always exists a family of 
closed internal systems {A(A),^}a£H given by abstract boundary conditions 
connected with the function t(-), such that the i?-matrix {i?(A)}AeH and the 
scattering matrix {S'(A)}agr of {L, L } can be expressed with the help of the 
eigenvalues and eigenfunctions of A(X) for a.e. A S R, cf. Theorem 15.51 This 
representation requires in addition that the internal Hamiltonians A(X) satisfy 
A(X) < Aq, which is always true if Aq is the Friedrichs extension of A. More- 
over, our general representation results also indicate that even for small energy 
ranges it is rather unusual that the i?-matrix and the scattering matrix can 
be represented by the eigenvalues and eigenfunctions of a single A-independent 
internal Hamiltonian A. 

,2 

As an application again the second order differential operators A = —-£p+ v 

,2 

and T = — + V from above are investigated and particular attention is paid 
to the case where the potential V is a real constant. Then the family {A(A)}agr 
reduces to a single selfadjoint operator, namely, to the Schrodinger operator in 
L 2 ((x;, x r )) with Neumann boundary conditions. In general, however, this is 
not the case. Indeed, even in the simple case where V is constant on (— oo,xi) 
and (x r ,oo) but the constants are different, a A-dependent family of internal 
Hamiltonians is required for a certain energy interval to obtain a representation 
of the i?-matrix and the scattering matrix in terms of eigenfunctions, see Sec- 
tion [672TT] The condition A (A) < Aq is always satisfied if Aq is chosen to be the 
Schrodinger operator with Dirichlet boundary conditions. Finally, we note that 
it is not possible to represent the i?-matrix and the scattering matrix in terms of 
eigenfunctions of an internal Hamiltonian with Dirichlet boundary conditions. 

The paper is organized as follows. In Section [2] we briefly recall some basic 
facts on boundary triplets and associated Weyl functions corresponding to sym- 
metric operators in Hilbert spaces. It is the aim of the simple examples from 
semiconductor modeling in Section 12.31 to make the reader more familiar with 
this efficient tool in extension and spectral theory of symmetric and selfadjoint 
operators. Section [3] deals with semibounded extensions and representations of 
Weyl functions in terms of eigenfunctions of selfadjoint extensions of a given 
symmetric operator. In Section [4] we prove general representation theorems for 
the scattering matrix and the i?-matrix of a scattering system which consists of 
two selfadjoint extensions of the same symmetric operator. Section[5]is devoted 
to scattering theory in open quantum systems, and with the preparations from 
the previous sections we easily obtain the abovementioned representation of the 
i?-matrix and scattering matrix of {L, Lq} in terms of the eigenfunctions of an 
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energy dependent selfadjoint operator family. In the last section the general re- 
sults are applied to scattering systems consisting of orthogonal sums of regular 
and singular ordinary second order differential operators. 



2 Boundary triplets and Weyl functions 
2.1 Boundary triplets 

Let $) be a separable Hilbert space and let A be a densely defined closed sym- 
metric operator with equal deficiency indices n±(A) = dimker(A* =p i) < oo in 
ft. We use the concept of boundary triplets for the description of the closed 
extensions of A in Sj, see e.g. [TI71 [TT1 [T2l [TO] . 

Definition 2.1 Let A be a densely defined closed symmetric operator in Sj. A 
triplet II = {H,r ,ri} is called a boundary triplet for the adjoint operator A* 
if Ti is a Hilbert space and To,ri : dom (A*) — * Ti are linear mappings such 
that the abstract Green's identity, 

(A*f,g) - (f,A*g) = (TJ^g) - (T f,T l9 ), 

holds for all /, g G dom (A*) and the mapping r := (£° ) : dom (A*) -^TL®TL 
is surjective. 

We refer to [IT] and [H] for a detailed study of boundary triplets and recall 
only some important facts. First of all a boundary triplet II = {Ti,T ,Ti} for 
A* always exists since the deficiency indices n±(A) of A are assumed to be 
equal. In this case n±(A) = dim7i holds. We also note that a boundary triplet 
for A* is not unique. 

In order to describe the set of closed extensions A C A* of A with the help of 
a boundary triplet II = {Ti,T ,Ti} for A* we introduce the set C(Ti) of closed 
linear relations in Ti, that is, the set of closed linear subspaces otTi®Ti. If 6 
is a closed linear operator in Ti, then will be identified with its graph Q(Q), 

@ = G{&) = Uq^ :ftedom(e) 

Therefore, the set of closed linear operators in Ti is a subset of C(Ti). Note 
that £ C(Ti) is the graph of an operator if and only if the multivalued part 
mul(6) := {ti £ Ti : ) <E 8} is trivial. The resolvent set p(&) and the 
point, continuous and residual spectrum er p (0), er c (0) and ay(O) of a closed 
linear relation & are defined in a similar way as for closed linear operators, cf. 
[13] , Recall that the adjoint relation 0* 6 C(Ti) of a linear relation in Ti is 
defined as 

0* := | W : (ti, k) = (h, k') for all W G ©} (2.1) 

and is said to be symmetric (selfadjoint) if C 0* (resp. = 0*). We note 
that definition (|2.ip extends the usual definition of the adjoint operator. Let 
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now 9 be a selfadjoint relation in Ji and let P op be the orthogonal projection 
in H onto H op := (mul (9)) 1 - = dom(9). Then 

is a selfadjoint (possibly unbounded) operator in the Hilbert space TL op and 9 
can be written as the direct orthogonal sum of 9 op and a "pure" relation 9oo 
in the Hilbert space Hoo '■= (1 — Pop)^ = mul 9, 

9 = 9 op ® 9 00 , 9 °° : ={U/J ^emuiejG^Woo). (2.2) 

With a boundary triplet II = {H.,Tq,Ti} for A* one associates two selfadjoint 
extensions of A defined by 

A Q := A* f ker(r ) and Ai := A* t ker(I\). (2.3) 

A description of all proper (symmetric, selfadjoint) extensions of A is given in 
the next proposition. 

Proposition 2.2 Let A be a densely defined closed symmetric operator in $j 
with equal deficiency indices and let H = {7i,ro,ri} be a boundary triplet for 
A*. Then the mapping 

Q^A e :=A* r T^Q = A* \ {/ G dom (A*) : (T f, IV) T e 9} (2.4) 

establishes a bijective correspondence between the setCiTL) and the set of closed 
extensions Aq C A* of A. Furthermore 

(A e y = a & , 

holds for any 9 G C{TL). The extension Aq in (12. 4|) is symmetric (selfadjoint, 
dissipative, maximal dissipative) if and only if 9 is symmetric (selfadjoint, dis- 
sipative, maximal dissipative). 

It is worth to note that the selfadjoint operator Aq — A* \ ker(ro) in (|2.3| 
corresponds to the "pure" relation 9oo = {(^) : h G 7i}. Moreover, if 9 is an 
operator, then (|2.4p can also be written in the form 

Aq = A* \ ker(ri - 9r ), (2.5) 

so that, in particular A\ in (|2.3|) corresponds to 9 = € [H]. Here and in 
the following [H] stands for the space of bounded everywhere defined linear 
operators in Ti. We note that if the product 9ro in (|2.5|) is interpreted in the 
sense of relations, then (|2.5p is even true for parameters 9 with mul (9) ^ {0}. 

Later we shall often be concerned with closed simple symmetric operators. 
Recall that a closed symmetric operator A is said to be simple if there is no 
nontrivial subspace which reduces A to a selfadjoint operator. By [23] this is 
equivalent to 

S) = clospan{ker(A* - A) : A e C\M}, 

where clospan{-} denotes the closed linear span of a set. Note that a simple 
symmetric operator has no eigenvalues. 
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2.2 Weyl functions and resolvents of extensions 

Let again A be a densely defined closed symmetric operator in with equal 
deficiency indices. A point A € C is of regular type if ker(A — A) = {0} and the 
range ran (A — A) is closed. We denote the defect subspace of A at the points 
A G C of regular type by N\ = kcr(A* — A). The space of bounded everywhere 
defined linear operators mapping the Hilbert space H into will be denoted by 
\H,Sf\. The following definition was given in [lOl lllj. 

Definition 2.3 Let A be a densely defined closed symmetric operator in 9), let 
n = {W,r ,Li} be a boundary triplet for A* and let A = A* \ ker(L ). The 
operator-valued functions 7 (-) : p(Ao) — > [7i,io] and M(-) : p(Aq) — » [H] defined 

by 

7 (A) := (r rA^) _1 and M(A) := L l7 (A), A € p(A ), (2.6) 

are called the 7-field and the Weyl function, respectively, corresponding to the 
boundary triplet II. 

It follows from the identity dom(A*) = ker(ro)+A / A, A G p(Aq), where as 
above Ao = A* \ ker(Lo), that the 7-field 7 (-) in (|2.6p is well defined. It is easily 
seen that both 7 (-) and M(-) are holomorphic on p(Aq), and the relations 

7(A) = (1 + (A - fx)(A - A)- 1 ) 7 (m), A, h e p(A ), 

and 

M(A)-M( M )* = (A-7i)7M*7(A), A,/x e p(A ), (2.7) 

are valid (see [H]). The identity (|2.7p yields that M(-) is a Nevanlinna function, 
that is, M(-) is holomorphic on C\M, M (A) = M(A)* for all A € C\R and 
3m(M(A)) is a nonnegative operator for all A in the upper half plane C+ = 
{A G C : 3m (A) > 0}. Moreover, it follows from that G p(Qm(M(X))) 
holds for all A G C\R. 

The following well-known theorem shows how the spectral properties of the 
closed extensions Aq of A can be described with the help of the Weyl function, 

cf. [mill]. 

Theorem 2.4 Let A be a densely defined closed symmetric operator in Sj and 
let {TL, Lo, Li} be a boundary triplet for A* with ^-field 7 and Weyl function M . 
Let Aq = A* \ ker(Lo) and let Aq C A* be a closed extension corresponding to 
some & G C(Ti.) via (|2.4[) - (|2.5[) . Then a point A G p{Aq) belongs to the resolvent 
set p(Aq) if and only ifO G p(0 — M(A)) and the formula 

{A @ - A)- 1 = (A - A)- 1 + 7 (A)(9 - M(A))"S(A)* (2.8) 

holds for all A G p(Aq) n p(Aq). Moreover, A belongs to the point spectrum 
<j p (Aq), to the continuous spectrum o~ c (A®) or to the residual spectrum o~ r (Ao) 
if and only i/0 £ o~i(@ — M(A)), i = p,c, r, respectively. 
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2.3 Regular and singular Sturm-Liouville operators 

We are going to illustrate the notions of boundary triplets, Weyl functions and 
7-fields with some well-known simple examples. 

2.3.1 Finite intervals 

Let us first consider a Schrodinger operator on the bounded interval (xi, x r ) C R. 
The minimal operator A in Sj = L 2 ((xi, x r j) is defined by 

^■■=-\ix^)Tx^ + < x ^ 

r i,y ^w^(( Xl ,x r )) ) (2.9) 

dom (A) := ^ feSj : f(x l ) = f{x r ) = 

{ (£/') = (£/0 (*r) = 

where it is assumed that the effective mass m satisfies m > and m, ^ € 
L°°((a;;, x r )), and that w e L°°((x;, x r )) is a real function, ft is well known that 
A is a densely defined closed simple symmetric operator in 9) with deficiency 
indices n+(A) = n-(A) = 2. The adjoint operator A* is given by 

(A*f)(x) = -~-L-^-f( x ) + v(x)f(x), 
K JJ{ ' 2dxm(x)dx Jy ' y > JK h 

dom (A* ) = {/ e Si : /, £/' G W 1 ' 2 ^, z r ))} . 
It is straightforward to verify that IT^ = {C 2 ,r ,ri}, where 

T0f = { f(*r) ) md Fl/ := 2 I - ft/') (* r ) J ' 

/ e dom (A*), is a boundary triplet for A* . Note, that the selfadjoint extension 
Aq := A* \ ker(ro) corresponds to Dirichlct boundary conditions, that is, 

dom(Ao) = {/€«:/, £/' € W 1 ' 2 ^, z r )), /(a:,) = /(av) = 0} . (2.10) 

The selfadjoint extension A\ corresponds to Neumann boundary conditions, i.e., 

Let (p\ and t/>a, A e C, be the fundamental solutions of the homogeneous differ- 
ential equation — \-^-^j^ u + vu = ^ u satisfying the boundary conditions 

<px{xi) = l, (^' x )(xi)=0 and 1> x ( Xl )=0, (^' x )(xi) = 1. 

Note that <p\ and V'a belong to L 2 ((xi,x r )) since (a;/,av) is a finite interval. A 
straightforward computation shows 

((A -\)-\f)(x)=vx(x) r Mt)f(t)dt+Mx) r <px(t)f(t)dt 

J Xi J X 

if\(x r ) 



lp\(x r ) 



1>x{x) / Mt)f(t)dt 



xi 
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for x G (xi,x r ), / G L 2 ((xi 7 x r )) and all A e p(A))- I n order to calculate the 
7-field and Weyl function corresponding to Ua = {C 2 ,r ,ri} note that every 
element fx G Af\ = ker(A* — A) admits the representation 

fx(x) = £,o<fx{x) + £,iip\(x), x£(xi,x r ), AgC, £ ,£iGC, 

where the coefficients £o>£i are uniquely determined. The relation 



To/a - 

yields 



1 \ 

^(av) ipx(x r )J Ui 



for ipx(x r ) 7^ (that is A G" <t(j4o)) and it follows that the 7-field is given by 
7 (A) :C 2 ^ L 2 ((x u x r )), 

(I) " ^k) ^"^^ ~ M-WxixrMo + M-)ti)- 

We remark that the adjoint operator admits the representation 



7 (A)*/ = 1 ( -C" (^ A ^) *l>x(x r) - j> \(y) fx(xr))f(y)dy\ 

f G L 2 ((x ; ,x r )). The Weyl function M(A) = Ti7(A), A G p(A ), then becomes 

1 (-<p\{x r ) 1 



M(A) - ^) V 1 -(^x)M 

All selfadjoint extension of A can now be described with the help of selfadjoint 
relations 9 = 9* in C 2 via (|2.4|) - ()2.5|) and their resolvents can be expressed in 
terms of the resolvent of Aq, the Weyl function M(-) and the 7-field 7(-), cf. 
Theorem 12.41 We leave the general case to the reader and note only that if 9 
is a selfadjoint matrix of the form 

e= (o °J' 



then 

dom (A @ ) = <J / = don.] (.4 

and 

(9-Af(A))- 1 = 



(^f')( X r) = -K r f(x r ) 



K r lpx(x r ) + {^^'x)( X r) 1 



V>A(x r )det(8-M(A)) \ 1 «#aOv) + <p\(x r ) 

9bviously the case ki — K r = leads to the Neumann operator A\. 
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2.3.2 Infinite intervals 

Next we consider a singular problem on the infinite interval (— oo, xi) in the 
Hilbert space Ri = L 2 {{— oo, x{)). The minimal operator is defined by 

(Tigi)(x) := —77-7- ttt-wW + U;0)sK a; ) I 

2 ax mi[x) ax 

where mi > 0, m>i,^- € 00, Xi)) and 17 e L°°((— 00, x/)) is real. Then 

XJ is a densely defined closed simple symmetric operator with deficiency indices 
n-(Ti) = n+(Ti) = 1, see e.g. [31] and [18] for the fact that Ti is simple. The 
adjoint operator T* is given by 

( T i*9i)(x) = -^"r— FTT^IX) + v i(x)9i(x), 
2 ax mi(x) ax 

dom(Tf) = {. 9; e «i : ffi, ^gl € ^((-oo.x,))} • 
One easily verifies that ilj; = {C, T , T^}, 

^off; := gi{xi) and T z l5i := - (-L-g[j ^ - 9i G dom ( T **)' 

is a boundary triplet for T* . Let <^a,; and ■i/'A,; be the fundamental solutions of 
the equation — ^-^j^j^u + v t u — Xu satisfying the boundary conditions 

<P\,l( x = (2^7^,/) OO = and ^\A x i) = °: (2^7^,/) ( x l) = !• 
Then there exists a scalar function m; such that for each A S C \ R the function 

a; 1 * ffA,«( x ) := <Pa,j( x ) - ttli(\)ipx,l( x ) 

belongs to L 2 ((— 00, xi)), cf. [39j. The function m; is usually called the 
Titchmarsh-Weyl function or Titchmarsh-Weyl coefficient and in our setting 
m; coincides with the Weyl function of the boundary triplet IIt, = {C, Tq, T^}, 
since 

r[g x , l = mi(\)r l g x ,i, gx,i € J\fx,i := ker(T* - A), AeC\K. 

An analogous example is the Schrodinger operator on the infinite interval 
(x r , 00) in Sif = L 2 ((x r , 00)) defined by 

(T r g r )(x) := -~ T^-j-gr(x) +v r (x)g r (x), 

2 ax m r (x) ax 

domfTV-L SSL- 9r^ r g' r eW^((x r ,oo)) 1 
dom(T r ).-| 9re ^. gAxr) ^ { j_ g , ){Xr) ^ ), 
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where m r > 0, m r , £ L°°((x r , oo)) and w r e L°°((a; r , oo)) is real. The 
adjoint operator T* is 

(T*g r ){x) = -~ l -—^-g r {x) + v r (x)g r {x), 

2 dx m r (x) dx 

dom(T;) = [g r £ : ffr, ^ € W 1,2 ((a; r , oo))} 
and n Tr = {C,T5,TJ}, 

T r g r := g r (x r ) and T\g r := ( — !— ) (x r ), g r £ dom(T r *), 

\2m r J 



is a boundary triplet for T*. Let ^A,r and ip\ tT be the fundamental solutions of 

1 d 1 d 

2 da; m r dx 



the equation ~\4z^-4zu + v r u — \u satisfying the boundary conditions 



cp\Axr) = 1, (a£;¥>A,r) (^) = aild ^\r{x r ) = 0, (^V^r) ( X r) = 1- 

Then there exists a scalar function m r such that for each A € C\R the function 

x >-> 5A : r(a;) := ¥>A,r(z) + m r (A)-0 A ,r(a;) 

belongs to L 2 ((x r ,oo)). As above m r coincides with the Weyl function of the 
boundary triplet IIt,. := {C, Tq, T^}. 

For our purposes it is useful to consider the direct sum of the two operators 
Ti and T r . To this end we introduce the Hilbert space 

£ := L 2 ((-oo, xi) U (x r , oo))=&i £ r . 

An element g £ .ft will be written in the form g = gi ®g r , where gi £ L 2 ((—oo,xi)) 
and g r £ L 2 ({x ri co)). The operator T = T\ T r in .ft is defined by 



— ^gi(x)+v l9l {x) 

-lE^E^W+MrW 



dom (T) = dom (T z ) © dom (T r ), 

and T is a densely defined closed simple symmetric operator in .ft with deficiency 
indices n + (T) = n_(T) = 2. The adjoint operator T* is given by 

{T9){X) -{ -\±^±9r{x) + v r9r {x) )> 

dom (T*) = dom (T,*) dom (T r *). 

One easily checks that n T = {C 2 ,T ,Ti}, T := (T[,,TS) T , Ti := {?[,T r 1 ) T , 
that is, 

To 3 =f ) and T l3 =if T^M^ 
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g e dom (T*), is a boundary triplet for T*. Note that T = T* \ ker(T ) is the 
restriction of T* to the domain 

dom (T ) = {g G dom (T*) : gi(xi) = g r (x r ) = 0}, 

that is, To corresponds to Dirichlct boundary conditions at xi and x r . The Weyl 
function r(-) corresponding to the boundary triplet TIt = {C 2 , To, Ti} is given 

A ~r W m _° A) ), A EpTO . 

3 Semibounded extensions and expansions in 
eigenf unct ions 

Let A be a densely defined closed symmetric operator in the separable Hilbert 
space $) and let {Tt,ro,Ti} be a boundary triplet for A* with 7-field and 
Weyl function M(-). Fix some 9 = 6* e C(H) and let A e C A* be the 
corresponding selfadjoint extension via (|2 .4[) . 

In the next proposition it will be assumed that Aq — A* \ ker(ro) and Aq 
(and hence also the symmetric operator A) are semi-bounded from below. Note 
that if A has finite defect it is sufficient for this to assume that A is semibounded, 
cf. Corollary O 

Proposition 3.1 Let A be a densely defined closed symmetric operator in j} and 
let {Ti, To, Til be a boundary triplet for A* with 7 -field 7(-) and Weyl function 
M(-). Let Aq be a selfadjoint extension of A corresponding to O = 0* £ C{Ti) 
and assume that Aq = A* \ ker(To) and Aq are semibounded from below. Then 
Aq < Aq holds if and only if 

ran( 7 (A)(8-M(A)) _1 ) Cdom(\/A -A) (3.1) 

is satisfied for all A < minjinf a(Ao), inf tr(Ae)} ■ 

Proof. Let Aq < Aq. From ([2~8| we get 

(A e - A)- 1 - (A - A)" 1 = 7 (A) (6 - M(A))"S(A)* > 

for A < minjinf <r(Ao), inf cr(Ae)} which yields 

(9- A/(A)) _1 > 0. 

By [TBI Corollary 7-2] there is a contraction Y acting from into Ti. such that 
(6 - M(A))~ 1/2 7 (A)* = Y(Aq - xr 1 ' 2 . 

Since A € K the adjoint has the form 

7(A) (6 - M(\)y 1/2 = (A e - X)-^Y*, 
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so that 

ran( 7 (A)(e-A/(A))~ 1/2 ) C dom(y/A e -X). 

Therefore 

ran( 7 (A)(e-M(A)) _1 ) C 

ran( 7 (A)(e-M(A)) _1/2 ) C dom(y/A e -\) 

and (|3.1[) is proved. 

Conversely, let us assume that condition (|3.1|) is satisfied. Then for each 
A < minjinf a(Ao), inf a(A&)} the operator 

FeW ■= v / ^^A7(A)(6-Af(A))" 1 (3.2) 

is well defined on Ti and closed, and hence bounded. Besides -Fq(A) we introduce 
the densely defined operator 

F e (A) = r (A e -A)- 1 / 2 , 
dom(F (A)) = {/ G £ : (A e - A)" 1 / 2 / G dom(A*)} 

for A < inf cr(ie)- 

It follows from IpTS]) . ^ = A* f ker(r ) and r 7(A) = 7 W that 

r (A e - A)" 1 = (6 - M(A)) _1 7 (A)* (3.4) 

holds for all A G p(A ) n p(A e ). Thus for A < min{inf cr(Ao), inf a(A&)} ([3~2jl 
becomes 

^e(A) = ^(rof^e - A)- 1 )* 
and together with (|3.3p we conclude 

F e (A) = r (A e - A)" 1 / 2 C (^/Z^A (r (A e - A)" 1 )*)* = (F£(A))\ 

This implies that -Fe(A) admits a bounded everywhere defined extension _Fe(A) 
for A < mimjmf a(A ),m£ a(A e )} such that F e (X)* = F e (X)* = F£(A). From 
(|3~4l) and M (A) = M(A)* we find 

r (r (A e - A)- 1 )* = (e - m(a))~\ x g p(a ) n P (A e ), 

so that for A < minjinf <j(Ao), inf cr(Ae)} 

(9 - M(A)) 1 = r (A e - A)- 1 / 2 ^^ (T (A e - A)- 1 )* 
= F e (A)F e (A)* >0. 

Using j2U) we find 

(Ae-A)- 1 > (Ao-X)- 1 
for A < minjinf <t(Ao), inf <r(Ae)} which yields Aq < Aq. □ 
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Corollary 3.2 Let A be a densely defined closed symmetric operator in Sj and 
let {7i,ro,ri} be a boundary triplet for A* with j-field and Weyl func- 
tion M(-). Assume that A has finite defect and that A = A* f ker(r ) is 
the Friedrichs extension. Then every self adjoint extension Aq of A in is 
semibounded from below and 

ran(7(A)(e-M(A)) _1 ) C dom(y/Ae - A) 

is satisfied for all A < minjinf cr(^4o)j m f cr(j4e)}- 

In the next proposition we obtain a representation of the function A i— > (0 — 
M(A)) -1 in terms of eigenvalues and eigenfunctions of Aq. This representation 
will play an important role in Section [5J 

Proposition 3.3 Let A be a densely defined closed symmetric operator in $j 
and let {7i,ro,ri} be a boundary triplet for A* with Weyl function M(-). Let 
Aq be a self adjoint extension of A corresponding to O = O* £ C(H) and assume 
that Aq — A* \ ker(To) and Aq are semibounded from below, Aq < Aq, and that 
the spectrum of Aq is discrete. Then the [H]-valued function A h- > (0 — M(A)) -1 
admits the representation 

oo 

(e-Af(A))- 1 = ^(A fc -A)- 1 (.,r ^)r ^ fe , \e P {A )np(A e ), (3.5) 
fc=i 

where {A^}, k = 1, 2, ... , are the eigenvalues of Aq in increasing order and {ipk} 
are the corresponding eigenfunctions. The convergence in (|3 . 5[) is understood in 
the strong sense. 

Proof. Let A < minjinf a(A ), inf o-(Aq)} and let E m , m € N, be the orthog- 
onal projection in onto the subspace spanned by the eigenfunctions {4>k}, 
k = 1, . . . ,m < oo, of Aq. Considerations similar as in the proof of Proposi- 
tion 13.11 show 

r ^ m 7(Ao)(e-A/(A )) _1 

= T Q (A e - A )- 1/2 S m VA e -Ao7(Ao)(e - M(A )) _1 
= FQ(X )E m FQ(X )* , 

where Fe(Ao) is defined as in (|3.3[) and Fe(Ao) S [Sj,H] denotes the closure. 
Hence we have 

lim r £ m7 (Ao)(e- M(A )) _1 = Fe(Ao)Fe(A )* = (0 - AT (Ao))" 1 

m— >oo 

in the strong topology. For A € p{A Q ) n p(A&) we conclude from the represen- 
tations 

(e-M(X))- 1 = v (T (A @ -x)- i y 

= Fe(A )(A e - A )(A e - A)- 1 F e (A )* 



R-MATRIX 



15 



and 

r £m7(A)(e - M(A)) 1 = Fe(Ao)(A e - A )(A e - A)- 1 S m F e (A )* 

that 

Hm r ^ m7 (A)(e - ai(x)) 1 = (e - a/(a)) _1 

m— *oo 

in the strong sense for all A £ p(A ) D p(Ae). 

Further, since the resolvent of admits the representation 

oo 

(A @ - A)- 1 = ]T(A fc - A) _1 (-, A e p(A e ), 

fc=i 

where the convergence is in the strong sense, we find 

m 

r (A e - \y 1 E m = J2(*k - A)- x (-, Vfc)roV-fc- 
*i=i 

For A € p(^4o) H /o(Ae) the adjoint operator is given by 
E m (T (A e - A)" 1 )* = E m ((e - Af(A)) _1 7 (A)*)* = E ml (\)(e - M(A))" 1 

m 

= ^(Afc-A)- 1 (- i r V fe )Vfe- 
fe=l 

Here we have again used (|2.8|l . A = A* f ker(r ) and r 7(A) = Replacing 
A by A and applying To we obtain from the above formula the representation 

m 

r 0J B m7 (A)(e-M(A)) _1 = ^(A fe -a)-^-,^)^ 

fc=l 

for all A € p(A ) n p(j4e). By the above arguments the left hand side converges 
in the strong sense to (0 — M(A)) -1 . Therefore we obtain (|3.5p . □ 

The special case 6 = £ [7Y] will be of particular interest in our further 
investigations. In this situation Proposition 13.31 reads as follows. 

Corollary 3.4 Let A be a densely defined closed symmetric operator in 9j and 
let {7i, Tq, r%} be a boundary triplet for A* with Weyl function M(-). Assume 
that A — A* \ ker(r ) and A\ — A* \ ker(ri) are semibounded from below, 
A\ < A , and that o~(A\) is discrete. Then the [H]-valued function A i— ► Af(A) -1 
admits the representation 

oo 

m(a)- 1 = ^(A-A,)- 1 (-,r v^)r v^, A G p(A ) n p{A x ), (3.6) 

k=l 

where {A&}, k = 1,2,..., are the eigenvalues of A\ in increasing order, {ipk} 
are the corresponding eigenf unctions, and the convergence in (|3.6p is understood 
in the strong sense. 
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Proposition 13.31 and Corollary 13.41 might suggest that the Weyl function M 
can be represented as a convergent series involving the eigenvalues and eigen- 
functions of the selfadjoint operator A . The following proposition shows that 
this is not possible if Ao is chosen to be the Friedrichs extension. 

Proposition 3.5 Let A be a densely defined closed symmetric operator in f) 
with finite or infinite deficiency indices and let {H,To,Ti} be a boundary triplet 
for A*. Assume that Aq = A* \ ker(ro) and A\ — A* \ ker(ri) are semi- 
bounded, that Aq coincides with the Friedrichs extension of A and that o~(Aq) is 
discrete. Then the limit 

m 

lim V(A-/x fe ) _1 (-, r i0fc) r i<fe: \ep(A ), 

m — >oo c — ' 

fc=l 

where {pk}, k = 1,2, ... , are the eigenvalues of Aq in increasing order and 
are the corresponding eigen] 'unctions, does not exist. 

Proof. We set 

Q(X):=T 1 (A -X)- 1 , Xep(Ao), (3.7) 

and 

G(A) : = r x <2(A)* = r^r^Ao - A)- 1 )*, a e P (A ). 

Taking into account the relation 

{At - A)- 1 = (A Q - A)" 1 - 7 (A)Af(A)~ 1 7 (A)*, A e p(A ) n p(M 

and (|2.6p we find 

g(A)= 7 (A)* and G(A) = M(A) A S p{A ) n p(Ai). 

Let m £ N, let E m be the projection onto the subspace spanned by the eigen- 
functions {4>k\, k — 1, . . . , m, and define 

Q m {\) := Q(X)E m and G m (X) :=TiE m Q(X)* , X e p(A ). 

With the help of 

OO 

(Ao-X)- 1 =^(/x fe -A)- 1 0,^ fe 
fe=i 

and (|3.7p we find the representation 

m 

G m (A) = ^( Mfc -A)- 1 (-,r 1 fc )r 1( /) fc , xe P (Ao)n P (A 1 ), 
fe=i 

and on the other hand 

G m (X) = Q m (X){A - X)E m Q(xy = 7 (A)*(A - A)£ m7 (A) 
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for Xe p(A )Dp(A 1 ). 

Let A G R, A < mm{mi a(Ao),mi a(Ai)}, and assume that there is an 
element rj G TL such that the limit 



lim G m (\)r]= lim V^(/ijfc — A) -1 (77, Ti<j) k ) T 1 cj} k 



(3.8) 



fc=i 



exists. Since for /i := 7(A)?] G A/a = ker(A* - A) 



(G m (A)*;, ry) = ((A - X)E ml (X)r,, 7 (A)r?) = \/A - A£ m /i 



we obtain from (|3.8|) that the limit lim„ woo HVA) — Ai? m /i|| exists and is finite. 
Therefore there is a subsequence {m„}, n € N, such that 



g := w- lim \/Ao — A E mn h and lim E mn h = h. 



Hence we conclude h G dom (V Ao — A) and g = \/Aq — Xh. But according to 
[IJ Lemma 2.1] we have dom (V-^o — A) DA/a = {0}, so that h — and therefore 
77 = 0. □ 



4 Scattering theory and representation of S and 
/^-matrices 

Let A be a densely defined closed simple symmetric operator in the separable 
Hilbert space Sj and assume that the deficiency indices of A coincide and are 
finite, n + (A) — n-(A) < 00. Let {TL,ro,Ti} be a boundary triplet for A*, 
Ao = A* \ ker(Lo), and let Aq be a selfadjoint extension of A which corresponds 
to a selfadjoint relation G C(H). Note that dim7i = n±(A) is finite. Let 
P op be the orthogonal projection in TL onto the subspace TL op := dom(O) and 
decompose 6 as in (I2.2[) . — Q op © 600 with respect to TL op © TLoo- The 
Weyl function M(-) corresponding to {TL, T , Ti} is a matrix-valued Nevanlinna 
function and the same holds for 

iVe(A) := (6 - Af(A)) -1 = (e op - Af p(A)) _1 P op , A G C\R, (4.1) 

where M op (A) = P op M(X)P op , cf. [25l page 137]. We will in general not distin- 
guish between the orthogonal projection onto TL op and the canonical embedding 
of H p into TL. By Fatous theorem (see [HI Hi]) the limits 

M(X + iO) := lim o M(A + ie) 

and 

N e (X + iO):= lim (6 — M(X + ie)) 1 
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from the upper half-plane exist for a.e. Ae K. We denote the set of real points 
where the limits exist by E M and E" 9 , respectively, and we agree to use a 
similar notation for arbitrary scalar and matrix-valued Nevanlinna functions. It 
is not difficult to see that 

iV e (A + »0) = (6 - M (A + iO))' 1 = (6 op - M op (A + iO))^ 1 P op , 

holds for all A e E A/ nS^ 8 and that M\(S M nS iVe ) has Lebesgue measure zero, 
cf. ]5l §2.3]. 

Since dim7i is finite by (|2.8[) 

dim (ran ((A - A)" 1 - (A - A)" 1 )) < oo, A <E p(A e ) n ^(Ao), 

and therefore the pair {Ao, ^4o} performs a so-called complete scattering system, 
that is, the u>aue operators 

W±(A e ,A ):= s- lini e^e^P^A)), 

t — >±oo 

exist and their ranges coincide with the absolutely continuous subspace Sj ac (Ae) 
of A e , cf. [HUGHES]- P ac (A Q ) denotes the orthogonal projection onto the 
absolutely continuous subspace fi ac (Ao) of Aq. The scattering operator Sq of 
the scattering system {Aq,Aq} is then defined by 

S e :=W + (A @ ,A )*W-(A @ ,A ). 

If we regard the scattering operator as an operator in $j ac (Ao), then Sq is 
unitary, commutes with the absolutely continuous part 

A a c :=A \ dom(A )n Sj ac (A ) 

of Aq and it follows that Se is unitarily equivalent to a multiplication operator 
induced by a family {Se(A)} of unitary operators in a spectral representation 
of Aq c , see e.g. [H Proposition 9.57]. This family is called the scattering matrix 
of the scattering system {Aq, A n }. 

In [1] a representation theorem for the scattering matrix {Sq(A)} in terms 
of the Weyl function M(-) was proved, which is of similar type as Theorem 14. II 
below. We will make use of the notation 

H MW := ran(9fm(M(A))), A e S M , (4.2) 

and we will usually regard TIm(\) as a subspace of Ti. The orthogonal projec- 
tion onto Wm(a) wu l be denoted by Pm(X)- Note that for A G p(A ) n R the 
Hilbert space 7~Cm(X) 1S trivial by (|2.7p . The family {Pm(a)}ag£ m °f orthogonal 
projections in Ti onto 'Hm(X)) ^ € is measurable and defines an orthogonal 
projection in the Hilbert space L 2 (M, d\,H). The range of this projection is 
denoted by i 2 (R, dX, H M {\))- Let P op and M op (A) = P op M(A)P op , A € S M , 
be as above. For each A 6 S M the space Hm(\) will also be written as the 
orthogonal sum of 

^M op (A) = ran(3fm(M op (A))) 
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and 

n M op (\) : = W M (A) © W Mop (A) = kcr(3m(M op (A))). 

The following theorem is a variant of [H Theorem 3.8]. The essential advan- 
tage here is, that the particular form of the scattering matrix {Sq (A)} imme- 
diately shows that the multivalued part of the selfadjoint parameter 9 has no 
influence on the scattering matrix. 

Theorem 4.1 Let A be a densely defined closed simple symmetric operator 
with equal finite deficiency indices in the separable Hilbert space and let II = 
{TI,Tq,Ti} be a boundary triplet for A* with corresponding Weyl function M '(•) . 
Furthermore, let Aq = A* \ ker(To) and let Aq be a selfadjoint extension of A 
which corresponds to = O pffi©oo G C(7i) via (|2.4p . Then the following holds. 

(i) The absolutely continuous part Aq c of Aq is unitarily equivalent to the 
multiplication operator with the free variable in L 2 (R, dX, Hm(\))- 

(ii) With respect to the decomposition TCm(\) = ^M op (A) © (A) ^ e sca ^~ 
tering matrix {Sq(X)} of the complete scattering system {Aq,Aq} in 
i 2 (R, dX, Hm(X)) ^ given by 

fSe ap (X) \ ± 

y M p (A) / 

where 

Se op (A) - In MopW + 2i v /Sm(M op (A))(e op - M^X))' 1 yjQm(M op (\)) 

and X G S M n T, N& , M op (A) := Af op (A + iO). 

Proof. Assertion (i) was proved in [4j Theorem 3.8] and moreover it was shown 
that the scattering matrix {S'e(A)} of the complete scattering system {Aq, Aq} 
in L 2 (R, dX,H.M(\)) has the form 

S©(A) = In MW + 2iV^m(Af(A))(e - Af(A))"V»ni(^W) e [W M( A)] 
for all A e E M n S^. With the help of 1(0) this becomes 



5 e (A) = I« M(A) + 2i A /3fm(Af(A))Pop(e p - M op (A)) P op ^Sm(M(X)). 



From the polar decomposition of ^Sm (M(X))P op , X G S M , we obtain a family 
of isometric mappings V(A), A G S M , from Hm (a) onto ran (ySm (Af (A))P p) 
dchncd by 

V r (A) v /3m(M op (A))a; := ^m(M(X))P op x 
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and we extend V(A) to a family V(X) of unitary mappings in Hm(X) ■ Note that 
V(A) maps ker(y / 3m (A/ op (A))) isometrically onto ker(P op ^/Qrn (M(A))). It is 
not difficult to see that the scattering matrix 

S&(\) := V(X)*Se(X)V(X), A e S M n E^ 9 , 

with respect to the decomposition Hm(\) = Hm op (\) © qn is of the form 
as in assertion (ii). □ 

We point out that the scattering matrix {Sq (A)} of the complete scattering 
system {Aq, Aq} is defined for a.e. A 6 1 and that in Theorem l4.1f ii) a special 
representative of the corresponding equivalence class was chosen. We also note 
that the operator -^/^ni (M op (A)) is regarded as an operator in 7iAf op (A)- 

Next we introduce the R-matrix {Rq(X)} of the scattering system {Aq, Ao} 
in accordance with Blatt and Weiskopf [5], 

Re(X) := i{l nMW - S e (X)) (I Hm{X) + Sq(X))- 1 (4.3) 

for all A G S M nS JVe satisfying -1 G p(S&(X)). Since Sq{X) is unitary it follows 
that i?e(A) is a selfadjoint matrix. Note also that 

S @ (X) = (Un Mm - ife(A)) + i?e(A))^ 1 (4.4) 

holds for all real A where R@(X) is defined. 

The next theorem is of similar flavor as Theorem 14.11 We express the R- 
matrix of the scattering system {Ae, A)} in terms of the Weyl function M(-) 
and the selfadjoint parameter 6 € C(H). Again we make use of the special space 
decomposition which shows that the "pure" relation part 0oo has no influence 
on the i?-matrix. 

Theorem 4.2 Let A be a densely defined closed simple symmetric operator 
with equal finite deficiency indices in the separable Hilbert space and let II = 
{H.,To,Ti} be a boundary triplet for A* with corresponding Weyl function M '(•) . 
Furthermore, let Aq = A* \ ker(To) and let Aq be a selfadjoint extension of A 
corresponding to G C(TC). Then for all A G £ M H E* 6 with 

ker(e op -SRe(M op (A))) ={0} 

the R-matrix of {Aq, Aq} is given by 

ReW = ^3m(M op (A))(e op -5Re(M op (A)))"V 3m ( M op(A)) 

with respect to H M (x) = T~Im op {\) © m, where M op (X) = M op (A + iO). 



Proof. It follows immediately from the definition (|4.3[) and the representation 
of the scattering matrix in Theorem 14.11 (ii), that the i?-matrix of {Ae,A)} 
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is a diagonal block matrix with respect to the space decomposition Ti-MtX) — 
^M p(A) © (A) an( ^ that the restriction of i?e(^) to ^, is identically 
equal to zero. 

Moreover, for every A S £ n S^ 8 it follows from the representation of the 
scattering matrix that 

v /3m(M op (A))(/ Wflfop(A) +5e op (A)) 

= 2{l nw) +^m(M op (A))(6 op - Af op (A))- 1 } v /3m(Af op (A)) 

= 2(e op - 3?e(M op (A)))(e op - A/ op (A)) _ y9m(Af op (A)) 

holds. If A e S M n S^ 9 is such that 9 op - SJe(M op (A)) is invertible, then we 
obtain 

y/$tm (M op (A)) (ln MopW + S @op (A)) 

= ^(©op - M op (A))(9 op - 5ie (Mop ( A) ) ) " 1 ^3ta (M op ( A) ) , 

so that 

2 l (9 op - Af op (A))-y3m(M op (A))(/ WMop(A) + 5 0op (A))^ 1 

= i(6 op - 5Re(A/ op (A)))" 1 v /3m(A/ op (A)). 

Finally multiplication by — ^/Sm (Af p(A)) from the left gives 

(^m„ pW --^eop (A)) (In + Se op (A)) 

= -iy / 9m(M op (A))(6 op - 3?e (M op (A)))~y9m(M op (A)) 

so that assertion (i) follows immediately from the definition of the i?-matrix in 
dO l. □ 



5 Scattering in coupled systems 

Let $j and & be separable Hilbert spaces and let A and T be densely defined 
closed simple symmetric operators in fj and respectively. We assume that 
the deficiency indices of A and T coincide and are finite, 

n := n+(A) = n-(A) = n + (T) = n_(T) < oo. 

Then there exist boundary triplets {T-L^Tq,Ti} and {7i,To,Ti} for the adjoint 
operators A* and T* , respectively, with fixed selfadjoint extensions 

A a := A* \ ker(r ) and T := T* \ ker(T ) (5.1) 



R-MATRIX 



22 



in and respectively, and dimW = n. The Weyl functions of {H, T , Ti} and 
{H,T , Ti} will be denoted by M(-) and t(-), respectively. Besides the spaces 
7~(-M(\), X G S M , (see (|4.2|l ) we will make use of the finite dimensional spaces 

H T{X ) = ran(Qm(T(A + iO))), A <= S r , 

and 

K(m+t)(A) = ran (3m ((Af + t)(A + iO))), A € E M+T D (S M n S r ). 

In the following theorem we calculate the 5 and i?-matrix of a special scat- 
tering system {L, Lq} in © M. in terms of the Weyl functions M and r. The- 
orem 15.11 is in principle a consequence of Theorem 14.11 and Theorem 14.21 cf . 
O Theorem 4.5]. We note that the coupling procedure in the first part of the 
theorem is similar to the one in [9]. 

Theorem 5.1 Let A, {H^Tq^x}, M(-) and T , {H,T ,Ti} ; r(-) be as above. 
Then the following holds. 

(i) The pair {L, Lq}, where Lq := Aq © To and 

L = A* © T* \ j/ © g € dom (A* © T*) : ^ + Ttf = } ' (5 ' 2) 

forms a complete scattering system in the Hilbert space $j © R and Lq c is 
unitarily equivalent to the multiplication operator with the free variable in 
L 2 (R,d\,H M (x)®n TW ). 

(ii) With respect to the decomposition 

W(M+r)(A) ®ft(M+T)(\) ( 5 - 3 ) 

ofriu(X) © ^r(A) ^ e scattering matrix {5(A)} of {L, Lq} is given by 
~ _ /5(A) \ ± 

- /-Hi 6 [«(M+t)(A) © "(M+r)(A)J ) 

\ n (M + T)(\) j 

where 

5(A) = Ih {m+t)W 

- 2 (M(A) + r(A)) (M(A) + t(A))"V^(^W + t(A)) 
and A € S M n E r n S( M+T ) _1 } M (A) := M (A + iO), r(A) = r(A + iO). 

(iii) For all A e S* f n S r n E^+^r 1 wift kcr(3f?c (M(X) + r(A))) = {0} the 
R-matrix of {L,Lq} is given by 

R(X) = 

/-V^m(M(A) +r(A))(Ke(M(A) +r(A)))"V» m (M(A) +r(A)) o\ 

V o 0^ 

tuif/i respect to the decomposition (15. 3p . 
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Proof, (i) Let L := A®T, so that T is a densely defined closed simple symmetric 
operator in the Hilbert space f) Clearly, T has deficiency indices n± (T) = 
2n, and it is easy to see that {H,r ,Ti}, where 

Fo(/ © 9) := (^° o Q , f i(/ ® 9) := and H:=H®H, 

f € dom(A*), g € dom(T*), is a boundary triplet for the adjoint operator 
L* = A* © T* in F) © .ft. Together with the selfadjoint operators Ao and To from 
(|5.ip we obviously have 

T :=L* fker(fo) - A Q (BT . 
It is not difficult to verify that 

6 := { ( ( ^)t) ■x ) yen\eC(H®H) (5.4) 

is a selfadjoint relation in 7i and that the corresponding selfadjoint extension 
T* \ T^^e in f) ffi £ via (f2T4j) coincides with the operator T in (JO}, cf. 
[3 . Since T has finite deficiency indices, T is finite rank perturbation of Lq in 
resolvent sense (cf. Theorem [23] and Section|4]), and hence {T, To} is a complete 
scattering system in S] © Moreover, as the Weyl function M(-) of {H, Tq, T±} 
is given by 

M(X)=( M ( X) r( ° A) ), Aep(To), (5.5) 

it follows from Theorem 14.11 (i) that the absolutely continuous part To c of To 
is unitarily equivalent to the multiplication operator with the free variable in 
T 2 (R, dX, Hu {x) ) = T 2 (R, dX, H M{X ) ® H T{X) ). 

(ii)-(iii) Note that the operator part O op of the selfadjoint relation 9 in (15. 4|) is 
defined on 

H op := dom (8) ={(x, x) T : x e U) 

and that 9 op = £ [Hop], cf. (|2.2|) . Next we will calculate the [H op ]-valued 
function M op (-), and in order to avoid possible confusion we will distinguish 
between embeddings and projections here. The canonical embedding of Ti op 
into H © Ji is given by 

1 

top : Hop ^Ti@Ti, — 

and the adjoint t* p g [HffiH,H op ] is the orthogonal projection P op from Tt(BTt 
onto H op , P op (u (&v) = -j^(u + v). Then we obtain 

Mop (A) = T op M(AKp = i(M(A) + r(A)), A € p(T ), 
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from (|5.5p . Now the assertions (ii) and (iii) follow easily from Theorem 14.11 (ii) 
and Theorem 14. 21 respectively. □ 

The case that the operator Aq has discrete spectrum is of particular im- 
portance in several applications. In this situation Theorem 15.11 reduces to the 
following corollary. 

Corollary 5.2 Let the assumptions and {L, Lq} be as in Theorem \5.1\ and as- 
sume, in addition, that o~(Aq) is discrete. Then the following holds. 

(i) Lq° is unitarily equivalent to the multiplication operator with the free vari- 
able in L 2 (R,d\,H rW ) 

(ii) The scattering matrix {S(X)} of {L,Lq} in L 2 (R, d\,Tl T (x)) * s given by 

5(A) = I HtW - 2i^m(r(A))(M(A) + r(A)) ~ V&n (r(A)) 
for A G S M nE r nl]( M + r ^ 1 ; where M{\) := M(X + iO), t(A) = r(\ + iO). 

(iii) For all A e n S r n E^+^r 1 with ker(M(A) + S?e(r(A))) = {0} the 
R-matrix of {L,Lq} is given by 

R(X) = -y/Sm (r(A)) (M (A) + Re (r(A))) " 1 y/Sm (r(A)) 

Proof. The assumption ct(Aq) — a p (Ao) implies 3m (JW(A)) = {0} for all 
A e S M . Therefore 

W(m+t)(A) = W T (A) and TCm(x) — {0}, A e E M . 

and the statements follow immediately from Theorem 15. II □ 

From relation (|4.4j) we obtain the next corollary. We note that this statement 
can be formulated also for the case when a(Ao) is not discrete. However in our 
applications we will only make use of the more special variant below. 

Corollary 5.3 Let the assumptions be as in Corollary \5.'A Then for all A £E 
S M n E r nEl M+r ) _1 with ker(Af(A) + 3?e(r(A))) = {0} the scattering matrix 
{iS'(A)} of {L,Lq} admits the representation 

SW = (iln rw + V^rn (r(A)) (M(A) + Rc (r(A))) ^ (r(A))) 

(iIn TW - (r(A)) (M(A) + 5Re (r(A))) " V»m (r(A))) ~* 

and, if, in particular, 3?e(r(A)) = 7 then 

SW = (iI HrW + V3m (r(A))M(A)" V9ta (r(A))) 

(i/*^, - v /3m(r(A))A/(A)-V3n 1 (r(A))) _1 . 
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Our next objective is to express the scattering matrix of the scattering system 
{L, Lo} in terms of the eigenfunctions of a family of selfadjoint extensions of 
A. For this let again r(-) be the Weyl function of {H, T , Ti}, let \i € £ T , and 
let {7i,ro,ri} be a boundary triplet for A* as in the beginning of this section. 
Then 3?e (t(/z)) is a selfadjoint matrix in H and therefore the operator 

A-Uo (t(m)) = ^* r kcr(rx + 5Re (r( M ))r ) (5.6) 

is a selfadjoint extension of A in .f), cf. Proposition 12.21 Note that by 
Theorem 12.41 a point A G p(^4o) belongs to /9(A_sft c ( r (^))) if an d only if 
G p(M(X) + 3fer(/x)) holds. The following corollary is a reformulation of 
Proposition 13.31 in our particular situation. 

Corollary 5.4 Let A, {ft,r ,ri} 7 M(-) and T, {W, T , Ti}, r(-) 6e as a&ove 
and assume o~(Aq) = a p (Ao) and that A is semibounded from below. For each 
/i£ E T with A_s% e < t <ij,)) < Aq the function A i— > — (5fte (V(/i)) + M(A)) -1 admits 
the representation 

oo 

-(M(A)+3?e(r( M )))" 1 =^(A fe M-A)- 1 (.,r ^M)r V fc M, 

fc=i 

where {Xk[tA}> fc = 1,2,..., are i/ie eigenvalues of the selfadjoint extension 
j4_sr c ( t (/j)) * n increasing order and ipk[lA we the corresponding eigenfunctions. 

Setting /i = A in Corollary 15.41 and taking into account Corollary 15.21 and 
Corollary [53] we obtain the following representations of the i?-matrix and scat- 
tering matrix of {L, Lq}. 

Theorem 5.5 Let the assumptions be as in Corollary \5.4\ Then for all A € 
S M n E r n S( M+T ) _1 with ker(M(A) + 5Re(r(A))) = {0} and A_ Uc(Tix)) < A 

the R-matrix and the scattering matrix of {L, L } admit the representations 

oo 

R(X) = - A)- 1 ( v /3m(r(A))-,r #[A])V»m WA))r ^ fc [A] 

k=i 

and 
S(\) = 

OO 

{ iJ Hr W -£(A*[A] - A)- 1 ( v /Sm(r(A))-,r V fe [A])^m(r(A))r ^[A]) x 

fe=l 

OO _ - 

+ E( A ^ W ~ A ) _1 ( V^m (r(A)) •, ToVfc [A]) (r(A))r Vfc [A]) " , 
fc=i 

respectively, where {Afc[A]}, fc = 1,2,..., are the eigenvalues of the selfadjoint 
extension A_$t e r T t\\) in increasing order and ipk[X] are the corresponding eigen- 
functions. 
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If SJe(r(A)) = for some A e £ r , then the operator A_sro(t(A)) m (|5.6D 
coincides with the selfadjoint operator = A* f ker(ri). This yields the next 
corollary. 

Corollary 5.6 Let the assumptions be as in Corollary \5.4\ Then for all A G 
S M flS T n E^+^r 1 sfte(r(A)) = 0, ker(M(A)) = {0} and A x < A Q the 
R-matrix and the scattering matrix of {L, Lq} admit the representations 

oo 

R(X) = J2(*k - A)~ 1 (v / 3m (t(A))-, ToVit) (r(A))r Vfc 
fe=i 

and 

oo 

5(A) =(tI WrW - £(A* - A)" 1 ( (r(A))-, r ^ fc ) V^m (r(A))r ^) 
fc=i 

oo _^ 

(i/H rW + E( Afc - A ) _1 (v 73 ™ (r(A))-, r Vfc) y/^tm (r(A))r ^ fe ) " , 
fc=i 

respectively, where {Afc}, k = 1,2,..., are i/ie eigenvalues of the selfadjoint 
extension A\ in increasing order and ipk are the corresponding eigenf unctions. 

Remark 5.7 The assumption A\ < A$ in Corollary 15.61 above is necessary. 
Indeed, let us assume that ^4o < A\ and that A\ is the Friedrichs extension. 
Let us show that in this case the sum 

oo 

^(Afc-AJ-^-.To^Oro^fc (5-7) 
fc=i 

cannot converge, where {A&} and {V'fc} are the eigenvalues and eigenfunctions of 
At. For this consider the boundary triplet {H^Tq,^}, T' = Fi and T[ = -T . 
Obviously A' {) = A* \ ker(r^) = A u A[ = A* \ ker(T[) = A and A' is the 
Friedrichs extension. By Proposition 13.51 we obtain that the sum 

oo 

^(A-Afc)- 1 ^,^)^ 

k=l 

diverges, where {Afc} and {ipk} are the eigenvalues and eigenfunctions of A' Q — 
A\. Using = — r one gets that the sum (|5.7p diverges. 

6 Scattering systems of differential operators 

In this section we illustrate the general results from the previous sections for 
scattering systems which consist of regular and singular second order differential 
operators, see Section [231 
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6.1 Coupling of differential operators 

Let the symmetric operators A = — + v an d 

T-T T - ( 1 d 1 d A ( 1 d 1 d 

1 r \ 2 dx mi dx J \ 2 dx m r dx ' 

in Sj = L 2 ((xi,x r j) and & — L 2 ((— oo, x{)) © L 2 ((x r , oo)) and the boundary 
triplets n A = {C 2 ,r ,ri} and n T = {C 2 ,T ,Ti} be as in Section [233] and 
and Section [2. 3. 21 respectively. By Theorem l5.1f i) the operator 

L:=A*®T*\S.f®gedom(A*®T*): r°j + Y°p = o} ^ 
is a selfadjoint extension of L — A © T in © .ft. We can identify ,f) © ^ with 

£ 2 ((z ; , x r )) © L 2 ((-oo, x/)) © L 2 {{x ri oo)) S L 2 (M). 

The elements /ffig in f)©ft, f € ft, g = gi®g r G ft, will be written as f®gi®g r . 
Here the conditions T / = T g and Tif = -Tig, / G dom {A*), g G dom (T*), 
explicitely mean 

S/fo) = f( x l) and /( x r) = g r {x r ), 

and 

-/') fa) = f — 9l) (xi) and (1/) (x r ) = (—A (x r ). 
m J \mi J \m J \m r ) 

Hence the selfadjoint operator (|€>. 1[) has the form 

L(f ®9i® 9r) = 

f-h£^£f+vf \ 

-l£^-s9i + vigi 

V -k£^:£9r + v r g r ) 

and coincides with the usual Schrodinger operator 

Hid „ . f „ ,o,-. . 1 



2 dx m dx 

where 



v\ /a 2 (i):/, Wgw 



1.2/ 



and 



m(x), x G (a;;, x r ) 

m(x) := ^ rrii(x), x G (— oo,a;/) 

m r (a;), x G (x r , oo) 

u(ic), a; G (xi, av) 

v(x) := { v t (x), x G (-oo,^) 

w r (a;), cc G (x r , oo). 
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The selfadjoint operator To = To, where Ao — A* \ ker(To) and 
T = T* f kcr(To), is defined on 

dom (To) = 

f®gi®g r e dom (A*) © dom (Tf ) © dom (T r *) : ^Z^^ Q 
and can be identified with the selfadjoint Schrodinger operator 

6.2 5 and .R-matrix representation 

It is well known that all selfadjoint extensions of the differential operator A 
in L 2 ((xi, x r )) have discrete spectrum. Hence according to Theorem 15.11 and 
Corollary 15.21 the selfadjoint Schrodinger operators T and To form a complete 
scattering system {T, To} in T 2 (R) and the scattering matrix {5(A)} is given 

by 

5(A) = I Ht{x} ~ 2iy/Zm(T(\))(M(\) + r(A)) "V^WA)) (6.2) 

for A G E M flE T n £( M + r )~ . Here M(-) is the Weyl function corresponding to 
the boundary triplet = {C 2 ,ro,ri} and 

»«*W-(f MA))' 

is the Weyl function of = {C 2 ,To,Ti}, cf. Section T2.3.2I It follows from 
[20] that for A 6 IT with Sm(r(A)) ^ the maximal dissipative differential 
operator 

A_ T(X) =A* rker^+rCA)^), 

that is, 

{A_ T{x) f){x) = -i—r/W + «(*)/(*). 

2 ax m(x) ax 

dom(A_ r(A) ) = { f e L 2 ((si s a:r)) : (^/') fo) = -mi(A)/(a;,) 
I («r) = m,(A)/(^) 

has no real eigenvalues, i.e. R C p(.A_ T m), so that each A e £ M = p(Ao) n M 
necessarily belongs to the set £( M+r ) by Theorem 12.41 Therefore the repre- 
sentation flUJ is valid for all A e {t e E T : Sm (r(f)) 7^ 0} n p(A ). Moreover, 
for AeS r with 3m(r(A)) = we have 5(A) = {0}. 

It is well known that the symmetric operator A given by (|2.9[) is semi- 
bounded from below and that the extension Aq = A* \ ker(ro), cf. (|2.10|) . 
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is the Friedrichs extension of A. In particular, this yields Aq < Aq for any 
other selfadjoint extension Aq of A. 

The selfadjoint operator A_jj e / T M\), A e £ r = E m ' n E mr , is given by 

(iU W A))/)(*) = + w(a ° /(:c) ' 

dom(A_ ffie(T(A)) ) = { / g ^((a^av)) : (^/') (a:,) = -He(m,(A))/(a:,) 
I Sr/')( a; r) = «e(m P (A))/( ! c r ) 

and clearly a(A_s% e ( T (A))) is discrete and semi-bounded from below for all A € 

Taking into account Theorem l5.5l it follows that the i?-matrix of {L, Lq} has 
the form 



i?^(A) = ^(A fc [A]-A)- 1 



k=l 



v /3?m(m r (A)) - J ' V>fc[A](s r ) 



yam(mi(A)) V» fc [A](a!i) 
v /Qm(m r (A))?/;fe [A] (a; r ) 

for all A G E r n E M with the property ker(M(A) + 5Re(r(A))) = {0} and 
3m(r(A)) ^ 0. Here {Afc[A]}, k — 1,2,..., denote the eigenvalues of the 
selfadjoint operator A_^ c ( t (a)) m increasing order and ipk[X] are the corre- 
sponding eigenfunctions. Furthermore we have again used R C p(A_ r ( A )) if 
9m(r(A)) 7^ 0, and moreover, R(X) — {0} if 3m(r(A)) = 0. 

The scattering matrix {^(A)} of {L, Lq} can be represented in the form 



fe=i 



\ /Sm (m;(A)) ^ fc [A] (xj ) 
v /Q?m(m r (A))?/;fe[A](x r ) 



fc=i 



VSm(m ; (A)) -\ /VfcW(a!i) 
v/3m(m r (A))7 ' \>[A](x r ) 



V'Sm (m ; (A))-0fe [A] (xj ) 



v ^/Sm(m r (A))^fc [A] (x r ) 
for all A e E r n E M with ker(M(A) + SJe (r(A))) = {0} and 3m (r(A)) ^ 0. 



6.2.1 Constant potentials vi and u r 

Let us assume that the potentials «;(•) and w r (") a s well as the mass functions 
m;(') and m r (-) are constant, that is, Vi(x) — vi G R, mi{x) — mi > for 
a; G (—00,0;;) and u r (x) = w r G K, m r (x) — m r > for a; G (u r ,oo). The 
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Titchmarsh-Weyl functions trij(-) and m r (-) can be calculated explicitly in this 
simple case, see [3]. One gets 



m/(A) = i\ — — and m r (A) = i\ 



2m; V 2m r 

for A G C+, where the square root is defined on C with a cut along [0, oo) and 
fixed by 3m (VA) > for A ^ [0, oo) and by \/A > for A G [0, oo). It is clear 
that 

E T = s mi n £ mr = E 
and it is not difficult to check 

{A G S r : 3m(r(A)) ^ 0} = (min-fa, v r }, oo). 

Furthermore 

Ke (mi (A)) 



2m; 1 — M 



0, A > vi, 



and 

SJe(m r (A)) = 



,-A 



A < u r , 

0, A > V r . 



If A € (max{u;, v r }, oo), then 5Re(r(A)) = and it follows from Corollary 
and the above considerations that the i?-matrix of {L,Lq} has the form 



y /Qm(mi(X)) ipk (xi ) 
v /3m(m r (A))V'fe (av ) 



(6.3) 



for all A G S M with the property ker(M(A)) = {0}. Here {A fc }, k = 1,2, . . . , 
denote the eigenvalues of the selfadjoint operator ^! in increasing order and ^fc 
are the corresponding eigenfunctions. Note that Ax is the usual Schrodinger 
operator in L 2 ((xi, x r )) which corresponds to Neumann boundary conditions, 
cf. (j2~TTj) , and that A G £ M has the property ker(Af(A)) = {0} if and only if 
A G p(A ) n p(Ax), cf. Theorem HH 
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Analogously the scattering matrix {5(A)} of {L, Lq} has the form 
5(A) =\ U H 



y /8m (m;(A)) V)fc (x; ) 
V3m (m r (A))V'fc(a; r ) 



fc=l 



\ /Sm (m;(A)) ^ fc (x; ) 
x /9m(m r (A))V'/c(a; r ) 

for all A £ (max{wi, zv}, oo) n p(A ) fl p(yli). 

The situation is slightly more complicated if A € (imn{vi, v r }, max{w;, v r }). 
Assume e.g. vi > v r and let A £ (v r ,vi). In this case 3m(r(A)) ^ 0, but the 
condition die (r(A)) = is not satisfied since 



5fte (m ; (A)) = - J— — - and 5Re (m r (A)) = 0. 
V 2m; 

The operator ^4-src(t(A)) is given by 

dom(A_ fic(T(A)) ) = I f G L 2 ((^,av)) : (^/') (*,) = yp£±f{xi) 

{ («r) = 

Since 

/0 \ 

v/9ta(r(A)) = ^ /a^^J, A€(«r,«l), 

the representations of the R and 5-matrix of {L, Lo} from the previous subsec- 
tions become 

oo 

R(X) = ^(A fc [A] - A)- 1 (V3m(m r (A))-,Vfc[A](x r ))73m(m r (A))^[A]( : r r ) 
fe=i 

and 



5 = j - E^i(AfcM - A)- 1 (yam (m r (A)) ■, MM(xr)) (m r (A))^ fc [A](3; r ) 
i + Er=iW A ]- A ) _1 ( \/9m (m r (A)) • , ^ k [A] (av)) y/Shn (m r (A))^ fc [A] (av) ' 



respectively, for A G (v r ,vi) n p(Ao) n p(-A-sR e (r(A))); see Theorem 15.51 Here 
{Afe[A]}, k — 1, 2, . . . , are the eigenvalues of the selfadjoint extension A^ c i T (X)\ 
in increasing order and ipk[^\ are the corresponding eigcnfunctions. 
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Remark 6.1 One might guess that the sum 




in the representation of the scattering matrix in (|6.3[) . where {Xk} and {ipk} are 
the eigenvalues and eigenfunctions of the Schrodinger operator with Neumann 
boundary conditions, can be replaced by the sum 

jJTj V V \2m ( Pk)\ X r)J J \ \2m C Pk)\ X r) J 

where {nk} and {4>k\ are the eigenvalues and eigenfunctions of the Schrodinger 
operator with Dirichlet boundary conditions. However, this is not possible since 
by Proposition 13 . 51 the last sum does not converge. We note that this can easily 
be verified by hand for the case v(x) = and m(x) = constant. 
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